Entanglement of two-mode Bose-Einstein condensates 
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We investigate the entanglement characteristics of two general bimodal Bose-Einstein condensates 
- a pair of tunnel-coupled Bose-Einstein condensates and the atom-molecule Bose-Einstein conden- 
sate. We argue that the entanglement is only physically meaningful if the system is viewed as a 
bipartite system, where the subsystems are the two modes. The indistinguishibility of the particles 
in the condensate means that the atomic constituents are physically inaccessible and thus the degree 
of entanglement between individual particles, unlike the entanglement between the modes, is not 
experimentally relevant so long as the particles remain in the condensed state. We calculate the 
entanglement between the two modes for the exact ground state of the two bimodal condensates 
and consider the dynamics of the entanglement in the tunnel-coupled case. 



I. INTRODUCTION 

In recent times, entanglement has come to be regarded 
as a physical resource which can be utilized to perform 
numerous tasks in quantum computation |lj] . This means 
that the creation and manipulation of entangled states is 
of significant interest in quantum information and com- 
putation. On the other hand, the study of the entangle- 
ment characteristics of various condensed matter systems 
has been proposed to provide new insights into quantum 
many-particle systems || |], ||, [| 0, ||, ||J . 

One extensively studied condensed matter system ]T^ , 
[l3| , [Til , 15 is that of a pair of tunnel-coupled Bose- 
Einstein condensates (BEC's). In the simplest model sys- 
tem, bosons are restricted to occupy one of two modes, 
each of which is a BEC. A dynamical scheme for engineer- 
ing many-particle entanglement between the particles in 
such a system has been proposed by Micheli et al [ [16] . 
However, since there is presently no definitive measure 
for entanglement between three or more subsystems, the 
amount of entanglement in the output state is not ana- 
lyzed quantitatively. Instead, this scheme aims to create 
states of a canonical form, whose entanglement content 
is based upon its inseparability. 

As argued in pp| , entanglement is only meaningful for 
multipartite systems whose Hilbert space can be viewed 
as a tensor product of two or more subspaces correspond- 
ing to physical subsystems of the system. As always, 
what one regards as an entangled state is, to some ex- 
tent, a matter of how this decomposition of the system 
into subsystems is performed. One person's entangled 
state is not the same as another's if they identify the 
subsystems differently. Entanglement can be said to be 
relative to the system decomposition |ll| . 

In the case of a bimodal Bose-Einstein condensate in- 
vestigated in (ll| H^|, the entangled subsystems were 
identified as the individual atoms in the condensate, and 
a mathematical measure of the multipartite entangle- 
ment proposed. By this measure certain states of the 
condensate were shown to be entangled. We argue here 



that the decomposition of the system in reference |T(| 
into subsystems made up of individual bosons is not 
physically realizable, due to the indistinguishibility of the 
bosons within the condensates Jl7], [ij] . We argue that a 
more physically relevant interpretation is to decompose 
the system into a bipartite system of the two modes. 
This idea is analogous to that of van Enk |l9| concerning 
the entanglement of electromagnetic field modes, as op- 
posed to the photons themselves. Since entanglement in 
bipartite systems is well understood, the entanglement 
between the modes of the system can be simply calcu- 
lated. We demonstrate that while the states created in 
reference Jl6| are indeed entangled, the entanglement has 
a different character. 

As a further extension, we consider another two-mode 
system, the atom-molecule BEC @ |lj H H @ H, 
p(| p^ . This systems has attracted significant interest 
since the entangled state is comprised of two chemically 
distinct components. We show that the entanglement be- 
tween the atomic and molecular modes can be calculated 
analogously to the entanglement between the two modes 
of the tunnel- coupled BEC's. 



II. THE SYSTEMS 

The situation where a large number of interacting 
bosons are restricted to occupy the same two-dimensional 
single-particle Hilbert space is known as the Josephson 
effect. The Josephson effect can be described as either ex- 
ternal, where the two single particle states, or modes, are 
separated spatially, or internal in which the two modes 
differ by some internal quantum number. Both the inter- 
nal and external Josephson effects can be described by 
the canonical Hamiltonian 11121 
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where a^,d^ denote the single particle creation oper- 
ators in the two modes (A and B) respectively, and 
Na — a \ a A, Nb — a B a,B are the corresponding boson 
number operators. The parameter Ej is the single-atom 
tunneling amplitude, A^i is the difference in the chemi- 
cal potential between the wells and K corresponds to the 
atom-atom interaction. Here we only consider K > 0, 
corresponding to a repulsive interaction between atoms. 
The total particle number, Na + Nb is a conserved quan- 
tity and is set to the constant value N. If we add the 
constant term 



K ( - - \ 2 
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the first term in the Hamiltonian (Q) becomes 
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as we expect for repulsive s-wave scattering. 

This Hamiltonian ([[]) is in fact a two-site version of the 
Bose-Hubbard model which describes bosonic particles 
with repulsive interactions, hopping through a potential 
lattice [^8], ^9|. In the Bose-Hubbard model, instead of 
two modes, there is an infinite lattice of potential wells 
(or modes) with coherent single-atom tunneling between 
nearest neighbor modes. 

A similar condensed matter system where there is the 
coupling of two BEC modes is that of atom-molecule 
Bose-Einstein condensate. In such a situation there exists 
coherent coupling between atomic and molecular BEC's 
respectively, which constitute the two modes of the sys- 
tem. The simplest Hamiltonian, recently studied by 
Vardi et al. 1 30 , which describes the atom- molecule BEC 
takes the form 



Ham = 2 + 2~ 
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where a' and w denote the creation operators for the 
atomic and molecular modes, respectively. Q is a mea- 
sure of the strength of the matrix elements for creation 
and destruction of molecules and S is the molecular bind- 
ing energy in the absence of coupling. The total atom 
number N atm = fi a + 2fii,, where n a = a^a, hb = b^b, 
commutes with the Hamiltonian, so is a constant of the 
motion. Both Hamiltonians (Fj) and (|J) have recently 
been shown by Zhou et al. |3l], ]3^] to be exactly solvable 
in the context of the algebraic Bethe ansatz. 

III. MANY-PARTICLE ENTANGLEMENT 

For qubits - two-dimensional systems represented by 
the states |0) and |1) - the canonical maximally entangled 
state is the Einstein-Podolsky-Rosen-Bohm (EPR) pair, 
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also known as a Bell state, in reference to the inequalities 
established by Bell J3j| . The tripartite analogue of this 
state is the Greenberger-Horne-Zeilinger-Mermin (GHZ) 
state 
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while the corresponding m-partite state is given by 
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These states are also known as m-particle Cat (m-Cat) 
states, in honour of Schrodinger's cat. For systems with 
three or more subsystems, while most definitely entan- 
gled, we cannot say whether the Cat states are maximally 
entangled. Since there is no definitive measure for arbi- 
trary multipartite entanglement, there is no clear notion 
of the structure of the maximally entangled sates in such 
systems. 

The d-dimensional analog of the 2-dimensional qubit 
is referred to as the qudit. For qudits, represented by 
the set of states {\i}} where i = 0, . . . , d — 1, a Cat state 
would be 



_L(|oo) + |(d-i)(<i-i))). 



(8) 



By the standard measure of entanglement for bipartite 
systems (the entropy of entanglement which is discussed 
in section [TV) this is not the maximally entangled state. 
While state (|) is entangled, a maximally entangled state 
is of the form 



d-l 



(9) 
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While Cat states are the canonical maximally entangled 
states for systems consisting of two qubits, for higher 
dimensions and number of subsystems, the maximally 
entangled states correspond to uniform distributions over 
the tensor product basis. 

In fuH , to determine the structure of the canonical 
entangled states, the system described by the Hamilto- 
nian (|l]) was decomposed into N subsystems consisting 
of the individual bosons, each with an internal degree 
of freedom described by a two dimensional Hilbert space 
spanned by the two states \A) and \B). In this descrip- 
tion, the system is viewed as a collection of N single- 
qubit subsystems. These internal degrees of freedom can 
be used to define a two mode description just as the po- 
larization degree of freedom of the electromagnetic field 
defines individual modes. In this case the annihilation 
and creation operators, appearing in equation ([l]), refer 
to single particle states distinguished by an internal de- 
gree of freedom rather than spatially localized single par- 
ticle states discussed in this paper. However this does not 
change our point of view regarding the lack of physical 
significance of entanglement at the level of single atoms. 
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In |16|] it was argued that the maximally entangled state 
in this case is the N-CaA, state, which is a coherent super- 
position state of all particles in mode A and all particles 
in mode B, i.e. 
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While it cannot be said that this is the maximally entan- 
gled state, it does indeed have some entanglement. 

However there is a problem with this choice of subsys- 
tem partitioning. By the nature of Bose-Einstein conden- 
sation, bosons within a condensate are indistinguishable. 
At no point can one make a physical measurement of the 
state of an individual particle in the condensate. For en- 
tanglement to exist between two systems the individual 
systems have to be distinguishable. While it is easy to 
imagine quantum measurements sensitive to individual 
particles, such operations could not be realized in the 
laboratory |l7|. While one can first remove individual 
particles from the condensate in order to measure them, 
the resulting state of the condensate is thereby changed, 
and it is unclear how the results of such measurements 
would reveal the multi-atom entangled state of the con- 
densate prior to the removal of the measured particles. 
This implies that the decomposition into individual bo- 
son subsystems is not physically realizable and while one 
can still write the Hilbert space of the system as a tensor 
product of the Hilbert spaces of individual bosons, this 
is not an appropriate description for realizable measure- 
ments upon the condensate. In other words, the system 
of coupled BEC's is best viewed as a bipartite entangled 
system rather than as a collection of N single-particle 
subsystems. 

Of course there is nothing to stop us from calculat- 
ing the entanglement between indistinguishable particles 
according to some measure. However entanglement is 
a physical resource that enables useful tasks in quan- 
tum communication and computation. In all such tasks 
it is necessary that the entangled subsystems be distin- 
guishable at some point in the protocols. For systems 
described by the Josephson Hamiltonian (Q) it cannot be 
said that there is physically useful entanglement between 
each individual boson when they exist in condensate. 



IV. ENTANGLEMENT BETWEEN THE TWO 
MODES 

Since the individual bosons are not physically acces- 
sible, distinguishable subsystems of the pair of tunnel- 
coupled BEC's described by ([!]), we need to consider 
other possible decompositions into subsystems if we are 
to investigate entanglement characteristics in this sys- 
tem. While we cannot measure which mode of the 
coupled BEC's a specific particle is in, the occupation 
number of a given mode is a physical observable. The 
two modes, be they spatially separated, or differing in 



some internal quantum number, are clearly distinguish- 
able subsystems. We can thus view the pair of coupled 
BEC's as a bipartite system of the two modes. It is rela- 
tively simple to investigate the entanglement between the 
modes since there is a unique measure of entanglement 
for two-component systems. Since the modes are distin- 
guishable the entanglement between them is accessible 
and thus potentially useful for some quantum informa- 
tion or communication protocol. This had been demon- 
strated by Dunningham et al. iQ who have proposed a 
scheme for entanglement swapping involving two pairs of 
tunnel-coupled BEC's. This is used to concentrate the 
entanglement between two modes. 

In this interpretation, while the entanglement involves 
many particles it is actually between the modes of the 
system. To illustrate this, consider the situation where 
we have just one particle in the system. In this scenario 
the modes can have occupation numbers of zero and one, 
so are spanned by the states |0) and |1). Consider the 
state 



i(|0)|l> 
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Clearly, with respect to the partition into modes, this 
single-particle state is entangled which implies we have 
entanglement with only a single particle. Analogous 
single-particle entanglement has been generated optically 



and used in a quantum teleportation protocol |36 



The state of each mode is characterized by its occupa- 
tion number. Because N is constant, a general state of 
the system \ijj) can be written in term of the Fock states 
by 
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where c n are complex coefficients, i.e., n bosons in mode 
A implies there are N — n bosons in mode B. 

The standard measure of entanglement of pure states 
of bipartite systems is the entropy of entanglement, which 
is the von Neumann entropy of the reduced density op- 
erator of either of the subsystems g, The reduced 
density operator of a subsystem is found by tracing out 
the other subsystem via the partial trace. If p is the den- 
sity operator describing some state of a bipartite system, 
the reduced density operator for subsystem A is defined 
by 



PA = Tr B (p) 



(12) 



where Tr# is the partial trace over subsystem B. The 
entropy of entanglement is then given by 



E(p A ) = -Tr(p A log(p A )) 
= -^AfclogtAfc) 



(13) 
(14) 



where the logarithm is taken in base 2, and {Afc} are the 
set of eigenvalues of the reduced density operator, p A . 
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The value of E varies between 0, for separable product 
states, to a maximum of log d (where d is the dimension of 
the Hilbert space of the subsystem) for maximally entan- 
gled states corresponding to a completely mixed density 
operator. 

The entropy of entanglement can be calculated from 
the reduced density operator of either of the subsystems 
without loss of generality. This follows from the Schmidt 
decomposition of pure states, which demonstrates that 
the eigenvalues of the reduced density operators of the 
two subsystems are identical ( pg. 109 of 0). 

Schmidt decomposition: For any pure state \il>) 
of a bipartite composite system there exist orthonormal 
states \ia) for subsystem A and orthonormal states \ib) 
for subsystem B such that 



U \IB) 



(15) 



where \ are non-negative, real numbers known as 
Schmidt coefficients, satisfying = 1. It is easy 

to see from the Schmidt decomposition that the reduced 
density operators for the two subsystems are, respec- 



tively, p A = J2i x i lu) (u| 



and p B = 



IB) (ib\ 



which have identical eigenvalues. 

Using the Fock basis, from equation ([ll]) the density 
operator describing a general state of the system is given 
by 



p = \i>) (V'l = 



c m c* n \m) \N-m) (n\(N-n\. (16) 



Taking the partial trace with respect to mode B yields 
the reduced density operator for mode A, 



PA 



trs (p) 

N 



= J2 c m c* n \m)(n\(k\N-m)(N-n\k) 
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|n) (n|. 



(17) 







Differentiating with respect to x n gives 

dL , 1 

— =»-lo gXn - — = 

which implies 

Xfi 2^ ln 2 
for all n. From the normalization condition 
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So a state with maximum entanglement will have coeffi- 
cients 



(22) 



Vn + T 



where n is some phase angle. This corresponds to a com- 
pletely mixed density operator, as expected for a state 
with maximal entanglement. Thus we can express the 
canonical maximally entangled state, \MES) for the sys- 
tem described by the Josephson Hamiltonian (Q) as 



1 N 

\MES) = -=== £ |n) \N - n) 



(23) 



From equation (|18[), the maximal entanglement is 



En 



N 
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= -log 
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log 
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1 



N + l 
log {N + l) 



(24) 



As mentioned previously, this is what is expected for 
the maximum entanglement, since the dimension of the 
Hilbert space of the individual modes is N + 1 (see page 
510 of |). 



From expression (|l^) we can see that the reduced den- 
sity operator in this case is diagonal in the Fock basis and 



the eigenvalues are simply Xi 



Thus the entropy 



of entanglement between the two modes of the coupled 
BEC's is given by 



E(pi) = - 



N 



n=0 



: log|c„| 



(18) 



To determine the maximally entangled state, expres- 



sion (18) can be optimized with respect to x n — 



by imposing the normalization condition 
with a Lagrange multiplier, p i.e. we maximize 

N 

L = -}^(x n loga;„ - px n ) + p. 

n=0 



o i u m 



= 1 



(19) 



A. Entanglement of the Ground State 

As mentioned in Section ||, the systems consisting 
of a pair of tunnel- coupled BEC's is the simplest sys- 
tem described by the Bose-Hubbard model and corre- 
sponds to a lattice potential with just two sites. For the 
Bose-Hubbard model, in the limit of an infinite lattice, 
there is a quantum phase transition where the ground 
state changes from superfluid phase to the Mott insula- 
tor phase |p8[ . Such a transition from the Mott insulator 
to superfluid phase was recently experimentally observed 
by Greiner et al. |]37| . 

In the Mott insulator state, particles tend to be lo- 
calized at the individual lattice sites with no phase co- 
herence across the lattice, whereas in the the superfluid 
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state, each atom is spread over the entire lattice and there 
exists long-range phase coherence across the lattice. This 
transition from the Mott insulator to the superfluid state 
occurs as the ratio of the coupling between lattice sites to 
the interaction strength increases. As long-range coher- 
ences in quantum systems are intrinsically linked to en- 
tanglement, it is of interest to quantify the entanglement 
in the system in relation to this transition. Since there is 
no measure for the entanglement in systems consisting of 
three or more subsystems, the two mode system here is 
the only Bose- Hubbard model for which we can currently 
give a complete description of the entanglement. 

Making the two modes identical (by setting the bias, 
Afi to 0), the Hamiltonian (|l|) was diagonalised numer- 
ically for increasing values coupling to interaction ratio, 
, and the entanglement of the ground state was deter- 



mined via equation (|l7|) . 
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FIG. 2: The variation of the entropy of entanglement of the 
ground state for increasing coupling to atom-atom interaction 



ratio, ^for N = 100. 




FIG. 1: The variation of the entropy of entanglement of the 
ground state for differing particle number, N, and increasing 
coupling to atom-atom interaction ratio, 

Figure [l] shows the results of this analysis for different 
values of the total particle number, AT. Since £j is the 
tunneling parameter, the larger its value the stronger the 
interaction between the modes of the system. As such, it 
is intuitive that that for no coupling, the entanglement 
between the modes is zero and as the coupling increases, 
the entanglement between the modes in the ground state 
increases. 

The entanglement asymptotically approaches a con- 
stant value as the ratio 4j ™> 00 , which is illustrated 
more clearly in figure ^|, which shows the results for 
N = 100. Now we consider the two extreme parame- 
ter values. Firstly for £j = 0, the ground state will have 
an equal, fixed number atoms in each of the two modes 
and is therefore the localized state, 



loc) 



(25) 



Clearly this state has zero entanglement. For K = 0, the 
Hamiltonian consists of just the tunneling term 



£ 



a qb + a B as 



(26) 



and it is easy to show that the ground state of such a 
Hamiltonian for a single particle is the bonding state 



71 ( fl 



10} |0) 



(27) 



where |0) |0) is the vacuum state. So in the ground state 
for K = each individual particle is in the bonding state 
and the state of the system is the TV-particle analogue of 
the bonding state 



I^) = 7Pm(^+4)" |o> io) 



(28) 



This state is the two-site analogue of the superfluid phase, 
with each atom being spread over the two modes. From 
equation ([l8|), the corresponding entanglement for this 
state is 



N 1 



. , 1 (N 
lQ g( ^JV 



(29) 



Thus for zero coupling, the ground state is the localized 
state (|25| ) which has no entanglement. As soon as the 
coupling begins to increase, the entanglement between 
the modes increases rapidly, as the occupation number 
of the modes is no longer exact and fluctuations in the 
phase decrease. As the tunneling amplitude continues to 
increase, the entanglement asymptotically approaches a 
maximum value, given by expression (|29|), which corre- 
sponds to the bonding state. In this state, the occupation 



6 



number fluctuations are large, resulting in phase coher- 
ence between the modes which is characterized by the 
high entanglement. 

Figure [3] shows that the maximum entanglement in the 
ground state is much less than the maximal entanglement 
of the system, log(7V + f ). 
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same independent of the value for D (> 1). The initial 
state prepared for the dynamical scheme of [jl6| is the 
iV-particle bonding state (|28|). According to expression 
(|29|), the entanglement between the modes in this state 
is actually greater than the entanglement in the final 
ideal output state, given by expression (30) with D = N. 
From this observation, it seems that the dynamical pro- 
cess outlined in |l6[| destroy entanglement between the 
modes. To gain a better understanding of this dynami- 
cal process with respect to the modal decomposition, we 
consider how the entanglement between the modes varies 
during the evolution. 



C. Dynamics of Entanglement 
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FIG. 3: The ratio of the entanglement of the bonding state to 
the maximal entanglement, E(p+)/ E m « for increasing par- 
ticle number, N. 



In studying the dynamics of the system of tunnel- 
coupled BEC's, we express the the Hamiltonian (Q) in the 
pseudo-angular momentum representation, introduced in 
p3[ , used in In this representation we define the 

three angular momentum operators 



(31) 
(32) 
(33) 
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However it appears the ratio of the entanglement of 
the two states remains finite as N — > oo. 



which have the canonical commutation relations 
[JxjJy] = iJz (and cyclic permutations). The Casimir 
invariant is easily seen to be 



B. Comparison of Entanglement 

The dynamical schemes developed in Ref.|]lq| aim to 
create states of the form (in the Fock basis) 



|Cat(D)) 



V2 ( 

N - D 



N + D 



N - D 



N + D 



(30) 



where D — Na — Nb- It was argued that the amount 
of entanglement in the state (^0|) is characterized by the 
"distance", D. For D = the state is separable and 
thus not entangled, while for D — N, |Cat(iV)) is equiv- 
alent to (|I(i| ) so is maximally entangled (which in this 
decomposition is still an incorrect assumption). However 
since motivation for the definition of these states was 
the decomposition of the system into individual particle 
subsystems, these states will have different entanglement 
characteristics when analyzed with the BEC modes as 
the subsystems. From equation jf8| ) it is clear that for 
D = the state (30) is separable and thus unentangled. 
However, for D > I, we have E{\CaX{D)) (Cat(Z))|) = I 
- the amount of entanglement between the modes is the 



J = 



N fN 



(34) 



In this way the tunnel-coupled pair of BEC's system is 
analogous to an angular momentum model with total an- 
gular momentum j — y. The Hamiltonian (0) (with 
A/i = 0) can thus be rewritten as 



H J2 = xj; - nj* 



(35) 



where we neglect constant energy shifts, x — if an d 
fl = ej. This Hamiltonian ( B 5 [) is slightly different from 
that defined in jlti) (and Jl3|] ). In these references, the 
J x term, which corresponds to the tunneling term, was 
added rather than subtracted (we assume only positive 
parameter values). However, this does not change the 
eigenstates of the system, but does reverse their order in 
terms of energies i.e. the ground state in the addition case 
is the highest excited state in the subtraction case. This 
means that the results for the dynamics from ||l6f can 
still be applied to Hamiltonian ([35]) by using a different 
initial state, as we will discuss below. 

In the angular momentum representation, states can 
be expanded in terms of the J z eigenstates, \j, m) z , where 
—j<m< j. In this basis there is no indication of 
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the underlying subsystem structure of the system. It is 
interesting to note that in terms of the Fock basis we 
have 

\j, m) z = \N — 2m) \N + 2m) . 
This implies that for any state 



N 



N 

= 5>„|j, n~N/2) z 



(36) 
(37) 
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meaning that the entanglement between the modes can 
be calculated from the coefficients in the angular momen- 
tum basis. 

In reference [[l6) , a semi-classical model of Hamiltonian 
( |35| ) was used to determine the optimal parameter values 
and time scale to create states of the form of expression 
@ with D = N (|Cat(A0)), from the evolution of a 
given initial state. In the angular momentum represen- 
tation 

\ Ca t(N)) = -j=(\j,-j) z + \j,j) z ) (38) 

From this analysis it was argued that using the critical 
parameter ratio 



2VL 

XN 



(39) 



in the evolution of the initial state over time, t c — 
ln(8N)/xN (where time is in units of h), could create 
states of the form d38|). 

The initial state given in |]l6|| was the maximal weight 



state of J x . This corresponds to the bonding state, (|28[ ). 
For the Hamiltonian above, the similar results will be 
achieved using the minimal J x weight state which is 



\^0)) = \j,-j) x ^-^=^ A -J B ) |0)|0). (40) 
In other words 



mt c )) = 



iH, 



|V(0)> » |Cat(JV)) 



(41) 



Using this parameter ratio and the initial state, the 
dynamics were investigated by numerically integrating 
the Schrodinger equation in the eigenbasis of J z . 

Figure ^ shows a plot of the evolution of the exact 
J z distribution, P(m,t) = \ z (j,m\ip(t)) \ 2 . This concurs 
with the results of showing that at time t c the final 
state is approximately of the form of the Cat state given 
in expression (38). 



The corresponding entanglement between the modes 
over the evolution is shown in figure |5|. The state at time 
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FIG. 4: The evolution of the J z distribution for j — 10 
(N = 20) with 2Q,/xN = 1. Note that at time t c w 2.5376 
the distribution is peaked at the two extreme m values, corre- 
sponding to an approximate Cat state as shown in jl6| . Again, 
the time is in units of h. 
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FIG. 5: The evolution of the entanglement. The dashed line 
shows the maximal entanglement, while the circle indicates 
the entanglement at time t c (time in units of h.) 



t c is not an exact Cat state, and as such the entanglement 
at t c was found to be greater than the initial entangle- 
ment. This differs from our earlier observation that the 
entanglement in the initial state is destroyed. This would 
be true if the state at time t c was exactly |Cat(7V)) (|3^). 

However, from figure ^ it is easy to see that the maxi- 
mum entanglement is not reached at this critical time but 
somewhat earlier, and at t c the entanglement between the 
modes is in the region of a local minima. As we can see, 
in the early stages of the evolution, the entanglement be- 
tween the modes approaches the maximal entanglement 
of the system. Relating this to figure || it is easy to 
see that the maximum entanglement is approached as 
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the original peak in the probability distribution flattens 
over the evolution, approaching a completely even distri- 
bution, corresponding to the maximally entangled state 
(^H). As the evolution continues, the distribution begins 
to peak at the extremes and the entanglement decreases. 

In terms of modal entanglement, the dynamical scheme 
proposed in [jl6| can still be used to create a close to 
maximally entangled state, over a shorter time period 
than for the creation of the inaccessible many-particle 
entangled state. 

It is of interest to note that the critical parameter ratio 
( |39| ) is the same as that found by Milburn et al. (l3) in 
regards to a transition in the dynamics from self-trapping 
to delocalization. For an initial condensate localized in 
one mode, when > 1, the condensate distribution 
will remain localized within the mode as it evolves. For 
< 1, the evolution results in a delocalization of the 

x 

condensate distribution between the two modes. 



where the partial trace has been taken with respect to the 
atomic mode and M is defined as above for even and odd 
total atom number, N atm . Thus the entropy of entangle- 
ment between the atomic and molecular modes is given 
by equation ([^J), the same expression as for the tunnel- 
coupled BEC's, where N, the total particle number, is 
replaced by M as defined above i.e. 



M 



E(p b ) = -^K| 2 log|d„ 



(45) 



n=0 



Since the dimension of the subspace of the modes is M, 
the maximally entangled states, analogous to (|23|), are 



\MES e 



1 M 

^=Y\2n)\M- 



n) , 



(46) 



for N atm even, and 



V. THE ATOM-MOLECULE BEC 

The atom-molecule BEC described by Hamiltonian (^) 
is a similar system to that of the tunnel- coupled pair of 
BEC's. In neither system can we consider the individual 
particles (the individual atoms and molecules) as sepa- 
rate, distinguishable subsystems but both consist of the 
coherent coupling of two distinct BEC's. In the atom- 
molecule BEC, the two modes of the system do not dif- 
fer spatially or by some internal quantum number but 
are rather two chemically distinct components. Nonethe- 
less, the determination of the entanglement between the 
atomic and molecular modes is analogous to the calcula- 
tions above for the tunnel- coupled BEC's. 

As before, the state of each mode is characterized by 
it's occupation number, however in the case of the atom- 
molecule BEC the set of Fock states spanning the Hilbert 
space of the system depends upon whether the total num- 
ber of atoms, Natm, is even or odd. In the case of an even 
N a trm a general state, \x) of the system can be expanded 
as 



lx) 



M 

E 

n=0 



d n \2n) \M - n) 



(42) 



where M = N atm /2, while for N atm 
state \(j>) can be expressed as 



M 

E 

n=0 



d n \2n+l)\M 



odd, the general 



(43) 



where in this case, M — (N atm — 1) /2 and the {d n } are 
complex coefficients defining the state. In analogy with 
expression (|lj) for the reduced density operator for the 
tunnel-coupled BEC's, the reduced density operator for 
a general state of the atom-molecule BEC is given by 



M 



Pb 



Eki 2 i m 



(M 



(44) 



\MESo. 



M 



n=0 



2n + l) \M 



(47) 



for N atm odd and will have entanglement log (M + 1). 

Fo llowing the same numerical analysis as in section, 
IV A figure ^ shows the results for the variation in the 
entanglement of the ground state of the atom-molecule 
BEC for differing values of the ratio of the parameters, 
^ and total number of atoms, N atm . To relate the en- 




100 



FIG. 6: The entropy of entanglement of the ground state of 
the atom-molecule BEC for increasing values of the ratio j| 
and atom number, N a tm- 

tanglemcnt structure of the ground state shown in figure 
^ to some physical properties of the system we need to 
consider other properties of the ground state for increas- 
ing total atom number and parameter ratio ^. Zhou et 
al. [ |32"| considered the two zero temperature correlations 
(n a ), the average atomic occupation number and the co- 
herence correlator, = (a^d^'b + b^ad). Figure shows 
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the results for the average atomic occupation number and 
the coherence correlator for the same parameter ranges 
used for the entanglement calculations. 




100 



FIG. 7: The average atomic occupation number (a) and the 
coherence correlator (b) for the ground state of the atom- 
molecule BEC. Note that both the average atomic occupation 
number and the coherence correlator have been scaled using 
the total atom number, N. 

We should note that the results here using direct nu- 
merical diagonalization of the Hamiltonian concur with 
the results found by Zhou et al. fl32| utilizing the exact 
solution. 

From figure f?]a it can be seen that the maximum en- 
tanglement in the ground state occurs where the average 
atomic occupation is comparative to the average molec- 
ular occupation. As indicated in |3^], the threshold cou- 
pling for the formation of a predominantly molecular 
BEC is ^ w 1A\/Natm- In the limit of large N atm , the 
threshold for the molecular BEC is in fact a quantum 
phase transition. Figure || shows the comparative results 
for the average atomic occupation number, the coherence 
correlator and the entanglement for N atm = 100. That 
the entanglement is not maximal at the quantum critical 
point is quite different to the behavior of the transverse 
Ising model, studied in references || |3f|. The entangle- 
ment characteristics of the transverse Ising model are, of 
course, much more complicated, since it consists of many 
distinguishable subsystems. In Q it was conjectured 
that, in the sense of entanglement sharing - how much 
two-party entanglement can be distributed amongst a 
given number of parties - the ground state was maximally 
entangled at the critical point. At the critical point the 
ground state saturates the bounds of entanglement shar- 
ing. While the ground state entanglement is not maximal 
at the critical point, the state is still strongly entangled. 
This would make intuitive sense, given that the property 
responsible for the long-range correlations in quantum 
phase transitions is entanglement. 




5/(ON IK ) 



FIG. 8: The average atomic occupation number, coherence 
correlator and the entanglement for the ground state of the 
atom-molecule BEC, for N a tm = 100. All three properties 
have been scaled with respect to their maximum value so as 
to compare the characteristics of these properties. 



The plots of the results for the entanglement and the 
coherence correlator share a common structure, however 
the maximum values occur for different parameter val- 
ues. This could mean that there is possibly another cor- 
relation that is more closely related to the entanglement 
between the atomic and molecular modes. 



VI. CONCLUSION 

We have argued here that in a system consisting of 
a pair of tunnel- coupled BEC's, the individual bosons 
within the condensates cannot be viewed as distinguish- 
able subsystems. Subsequently entanglement in this sys- 
tem should not be viewed as between the individual 
bosons. A more physically relevant description of this 
system is as a bipartite system where the subsystems are 
the two modes. Using this description we have analyzed 
quantitatively the entanglement between the two modes 
in the ground state of the coupled BEC's and its relation 
to the Mott insulator to supcrfluid phase transition. This 
idea was extended to consider the entanglement between 
the atomic and molecular modes in an atom-molecule 
BEC. 

On top of this, we have demonstrated that the dynam- 
ical scheme of [|l6| , argued to be viable with current state 
of the art technology, can be used to create a highly en- 
tangled state between the modes of the BEC system, over 
a smaller time-scale. 

As mentioned earlier, the amount of entanglement de- 
pends upon how the system is decomposed. In section 
IV C, it was shown that the tunnel-coupled two- mode 



system can be viewed as a pseudo- angular momentum 
system - a single audit. In this description - viewing the 
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solely as a single qudit - it appears that there is no en- 
tanglement present in the system. Entanglement is only 
seen when the system is viewed in terms of subsystems, 
in this case, the two modes. In other words, the entangle- 
ment cannot be characterized when we neglect informa- 
tion about the underlying subsystems and only consider 
properties of the system as a whole. 

A possible way to create entanglement between 
individual bosons in the tunnel-coupled system would 
be to engineer some state within the condensate traps, 
then free the particles (see |3^, [lo] for examples of 
this procedure applied to other BEC systems). Once 
the bosons are free from the condensate they become 
distinct allowing entanglement to form between them. 



However, whilst the bosons remain in condensate they 
are indistinguishable and cannot become entangled with 
each other. 
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